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L^: ABSTRACT 

(N' 

Q^ ■ We have observed a region of nebulosity first identified as starUght scattered 

OO . by interstellar dust by Sandage (1976) using the GALEX ultraviolet imaging 

telescope. Apart from airglow and zodiacal emission, we have found a diffuse 
UV background of between 500 and 800 photons cm~^ sr~^ s~^ A~^ in both 
the GALEX FUV (1350 - 1750 A) and NUV (1750 - 2850 A) bands. Of this 

O ■ emission, up to 250 photons cm~^ sr~^ s~^ A~^ is due to H2 fluorescent emission 

in the FUV band. The remainder is consistent with scattering from interstellar 

S^ . dust with forward scattering grains of albedo about 0.4. These are the highest 

spatial resolution observations of the diffuse UV background to date and show 
an intrinsic scatter beyond that expected from instrumental noise alone. Further 
modeling is required to understand the nature of this scatter and its implications 
for the ISM. 

Subject headings: ultraviolet: ISM, ISM: dust 



INTRODUCTION 



Ever since the f irst o bservations of diffuse ultraviolet radiation by lHayakawa et al.l (119691 ) 



and iLillie fc WittI (119761 ) . there has been an effort to understand its distribution and its 
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origin. Unfortunately, because of the difficulty of the observations and the faintness of the 
background, many of the early observations were conspicuous more by their disagreements 
than by the light they shed on the topic. T he s tate of the ob servations and theories before 
1990 have been reviewed by iBowyerl (jl99ll ) and iHenryl (jl99ll ). 



There has been signific ant progress in more rece nt years, particularl y in th e far ultra- 
violet (A < 1200 A) where iMurthv et all Jl999h and iMurthv fc Sahnowl JiooJ have used 
spectroscopic data from the Voyager and FUSE {Far Ultraviolet Spectroscopic Explorer) 
spacecraft, respectively, to trace the radiation field over many different locations in the sky. 
There have also beei i a number of ob servations at longer wavelengths, most recently by the 
SPEAR instrument (JRyu et al.ll2008l . and references therein), but no systematic study of the 
UV background. The Galaxy Evolution Explorer ( GALEX) offers us the opportunity to ex- 
tend coverage of the diffuse background to a significant fraction of the sky with a sensitivity 
of better than 100 photons cm^^ sr^^ s^^ A^^. In this work, we will r eport on one such 
observation: that of the nebulosity observed near M82 by ISandagd (Il976l ). as a template for 
our further work with a much larger data set. 

This cloud is a t a high Galactic latitude (38°) with few nearby stars and was identified 
by ISandagd (119761 ) as a canonical high latitude dusty cloud illuminated by the Galactic 
interstellar radiation field (ISRF). Our expectation was that we would be able to differentiate 
between starlight scattered from the Galactic cloud and extragalactic light which would be 
shadowed by the foreground cloud. The observations presented here are the first to probe the 
diffuse UV background at a spatial resolution comparable to other surveys of dust emission, 
notably the JR. 



OBSERVATIONS 



A full descriptio i i of th e GALEX instrxiTaent and mission is given by lMartin et al.l (120051 ) 
and iMorrissey et al.l (120071 ) with those parameters relevant to our observation listed in Table 
[U Although we had originally intended for an exposure time of 15,000 seconds in each of 
the far ultraviolet (FUV) and near ultraviolet (NUV) bands, power problems in the FUV 
high voltage power supply and subsequent re-observations of the region gave us about 15,000 
seconds exposure time in the FUV and 35,000 seconds in the NUV. 

A single GALEX observation is actually made up of a number of visits, each of approx- 
imately 30 minutes in length (limited by the duration of orbital night at this altitude) and 
possibly separated by several months. Although GALEX is a photon counting instrument, 
the standard GALEX pipeline combines all observed events to produce a single image per 
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channel per visit. We coadded the data from each visit to produce a single view of the 
observed region in each of the two channels (Fig. [1]). 

These images contain both stars and the diffuse background but, because our focus is 
on the diffuse background, we have eliminated all the stars by cutting out a box around each 
source listed in the GALEX merged catalog and then averaging the data into bins of 2' x 
2' in size. Note that we eliminated the edge effects seen in all GALEX observations before 
binning. These data represent the total diffuse radiation from all sources in this direction 
and in this observation. It is difficult for any imaging experiment to separate the different 
components of the diffuse radiation and, in the remainder of this section, we will discuss the 
possible sources of the background. It is important, however, to note that the direct stellar 
contribution (the sum of the counts from all the stars in the image) to the total signal in 
each channel is less than 7%. 



2.1. Instrumental Effects 

Photon counting instruments such as the delay-line detectors used in GALEX are intrin- 
sically low noise and the instrumental dark count, due largely to cosmic rays and other fast 
particles, is only o n the order of 20 and 60 counts per second over the FUV and NUV detec- 



tors, respectively (iMorrissey et al.ll2007l ). corresponding to an effective diffuse background of 
about 5 photons cm~^ sr^^ s^^ A~^. There are other instrument related effects which may 
give rise to artifacts in the background, most important of which is the flat field already 
applied to the data. Specifically to mitigate possible variations in the fiat field, GALEX 
observations are dithered over a 90 ^^ spiral leaving an empirical uncertainty of about 5% in 



the fiat field (IMorrissey et al.l 120071 ). We would expect an even smaller variation in our data 



where we further bin over 2' square boxes. 



2.2. Airglow 

A significantly greater part of the emission is due to the resonantly scattered O I lines 
at 1304, 1356, and 2471 A from the Earth's ionosphere with a contribut ion of between 10 



and 200 photons cm"^ sr"^ s~^ A"^ in each of the two GALEX bands flBoffi et ah 1120071 ). 
Uniquely amongst the contributors to the diffuse radiation, the airglow will vary over the 
course of a single visit. Indeed, plotting the total count rate from the GALEX TEC counter 
(Fig. [2]) shows that the airglow in a given visit is a function of the local time, with a 
minimum at local midnight. However, there is also a variation between visits which appears 
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to be roughly correlated with the level of solar activitjo (Fig. [3]) as would be expected from 
the resonantly scattered Solar photons in the Earth's ionosphere. There are two points in the 
NUV channel which blatantly violate this general trend. These are the only two in which the 
Sun angle is less than 90° and it is likely that the increased baseline is due to instrumental 
scattering of off- axis Solar light. These were not observed in the FUV because of the power 
issues at the time. 

We have adopted an empirical method to estimate the airglow in each observation by 
setting the airglow contribution to at local midnight on Jan. 3, 2006. The total airglow 
contribution to each visit is then the integral under the curve. These values are tabulated 
for the FUV channel in Table [2] and imply an average contribution to the total diffuse 
background from the airglow of about 85 photons cm~^ sr"^ s^^ A^^. A similar exercise for 
the NUV channel yields an average contribution of 120 photons cm~^ sr^^ s^^ A^^. Given 
the uncertainty in the baseline, we estimate that the total airglow contribution in our data 
is between 100 and 200 photons cm~^ sr~^ s~^ A~^, consistent with the levels observed at 
the slightly lower Hubble Space Telescope altitude of 600 km. 



2.3. Zodiacal light 

Moving out from the Earth, the next major contributor to the diffuse radiation field 
is the zodiacal light, sunlight scattered by interplanetary dust grains. The visible light 
dis tribution of the zodi acal light has been mapped as a function of helioecliptic coordinates 



by iLeinert et al.l (119981 ) and we have estimated the contribution of the zodiacal light to the 
NUV band by assuming that the UV distribution is the same as that in the visible with a 
color close to 1; i.e., the ratio between the zodiacal light and the solar spectrum is the same 
at all wavelengths. With these assumptions, the effective contribution of the zodiacal light 
is less than 10 photons cm~^ sr~^ s~^ A~^ in the FUV, where the solar spectrum vanishes, 
and 400 - 470 photons cm~^ sr~^ s^^ A^^ in the NUV, depending on the date of the visit. 

The main u ncertainty in p r edicti ng the level of the zodiacal light is that the color may 



differ from unity. iMurthy et al.l (119901 ) found a dependence of the color on the ecliptic latitude 
varying from 0.6 on the ecliptic plane to 1.2 at an ecliptic latitude of 60°, and 1 at the 50° 
latitude of these observations. Taking this into account, we estimate a 10% uncertainty in 
the zodiacal light or about 50 photons cm^^ sr^^ s^^ A^^. It should be reemphasized that 
there is essentially no zodiacal light contribution to the FUV channel. 



^archived at http://www.cixlc.com by Jan Alvestad 
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2.4. Ho Emission 



This GALEX field is included in one of tlie regions observed by lMartin et al.l (19901) with 



the Shuttle-based Berkeley Extreme Ultraviolet/Far Ultraviolet Shuttle Telescope {BEST) — 
their Target 6, a scan over Ursa Major. They detected emission from the Werner bands of 
molecular hydrogen over the entire observed region and speculated that it was due to an H2 
halo extendin g well outside the Sandage cloud. Extrapolating into the GALEX bands with 



the models of iMcCandlisd ( 120031 ) . this emission implies a uniform brightness of 135 photons 
cm~^ sr~^ s~^ A~^ over the FUV field with no contribution in the NUV channel. We will 
discuss the H2 contribution to our observed emission below but note that we have found the 
level of the emission to vary between and 250 photons cm~^ sr"^ s"^ A^^ over the field. 



2.5. Total Emission 

The contribution of each of the sources of diffuse emission listed above to the total 
signal is tabulated in Table [3] implying a uniform contribution of about 90 photons cm~^ 
gj^-i g-i ^-1 jj^ ^Yie FUV channel and 570 photons cm~^ sr~^ s^^ A~^ in the NUV, with an 
uncertainty of up to 100 photons cm~^ sr~^ s~^ A~^ in the baseline. We have subtracted 
these uniform baselines from each channel and the resultant images are shown in Fig. |H 
The remaining signal includes H2 fluorescent emission, starlight scattered by the interstellar 
dust in the field and an extragalactic component, if any, and will vary across the image. We 
note that the foreground emission comprises about 15% of the total signal in the FUV and 
close to 50% of the total signal in the NUV, where the zodiacal light is of roughly the same 
magnitude as the Galactic sky background. 

The uncertainty in the baseline is about 15% of the total signal but it should be em- 
phasized that this is not a pixel by pixel variation but rather a global uncertainty in the 
magnitude of the signal. The relative pixel-to-pi xel uncertainty is large ly determined by the 



uncertainty in the fiat field and is less than 5% (JMorrissey et al.l 120071 ). 



3. RESULTS AND DISCUSSION 

The emission in the FUV channel is strongly correlated with that in the NUV channel 
(Fig. \5^, as would be expected from starlight scattered by interstellar dust. The illuminating 
ISRF is spectrally flat in this region and, purely coincidentally, the optical depth of normal 
Galactic dust is the same in both channels with the FUV rise in the extinction curve balanced 
by the 2175 A bump in the NUV. However, a closer glance at the figure shows that the ratio 
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between the two bands (FUV/NUV) increases with the FUV emission (Fig. [6]), suggesting a 
source which contributes to the FUV band but not the NUV. This is most probably due to 
molecular hydrogen fluorescence and we have estimated its magnitude and distribution by 
assuming that the ratio between the two bands to be fixed at 0.8 — the minimum observed 
ratio — and attributing the remainder to H2 emission (Fig. [7]). The absolute level may be 
uncertain by about 100 photons cm~^ sr^^ s^^ A^^ as discussed above, but this will not 
affect the morphology of the emission. 

The level of the H2 varies to a maximum of about 250 photons cm~^ sr^^ s~^ A^^ over 
the f ield with no correlat ion (Fig. [8]) with the IRAS 100 ^tm emission or, by exte nsion, with 
CO (jWeiland et al.lll986l ). These results are consistent with lMartin et al.l (Il990l ) who found 
from their UVX/BE ST results that the H2 emission extended much beyond the CO contours 
and suggested that it was due to an extended halo of molecular hydrogen. However, we do 
find a spatial variability in the emission to which they were not sensitive because of the speed 
of their scanning and the lower spatial resolution of the instrument. The remaining emission 
can be attributed to dust scattered radiation of 450 - 650 photons cm~^ sr~^ s~^ A~^ in the 
FUV and 550 - 750 photons cm~^ sr~^ s~^ A~^ in the NUV, with a spatial variation greater 
tha n that due to photo n noise alone. Twenty years after UVX, it is difficult to correlate 
the iMartin et al.l (Il990l ) data with ours but it is clear that the two sets of observations have 
yielded consistent results. 

Our UV emission is fia t relative to the 100 ^m emission (Fig. [9]) in contrast to the 
strong correlation found by iHaikala et al.l (119951 ) in FAUST observations of another cirrus 
cloud (G251.2+73.3). Both observations are of isolated clouds illuminated only by the ISRF 
but the optical depth is much greater in the Sandage region. Because the UV emission is 
due to scattering from the surface of the cloud while the IR is due to thermal emission from 
the entire volume of the cloud, because of the low optical depth in the infrared, we would 
expect a correlation between the UV and IR emission at low optical depths (in the UV) with 
the UV emission saturating as the 100 /im emission increases. This is essentially what is 
seen in Fig. [9l with the caveat that a more detailed modeling of both regions is required. 



We ha ve used our standard model for scattering from interstellar dust (ISujatha et al. 



20051 . I2OO7I ) to derive the optical constants of the interstellar grains in this direction. This 
procedure uses the position and spectral type of the stars from the Hipparcos catalog to pre- 
dict the radiation field at the location of the scattering dust, assumed to be at a distance of 
120 D C. The stellar radiation is convolved with a Henyey-Greenstein (JHenyey fc Greenstein 
19411 ) scattering function for the grains to obtain the distribution and magnitude of the 
dust-scattered light in this direction. In our previous stu dies of the Ophiuchus region 
(ISujatha et al.l l2005l ) and the Coalsack ( Sujatha et al.l 120071 ). the scattering was from an 
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optically thin region with a small number of illuminating stars. Here, many more stars con- 
tribute to the radiation field with an optically thick dust cloud and it was not practicable 
to run a Monte Carlo multiple scattering model for the entire region. Therefore we assumed 
single scattering throughout the region but applied a correction for multiple scattering based 
on limited Monte Carlo runs. 

Two representative model fits (A and B) to the NUV data with the airglow emission 
and zodiacal light subtracted are shown in Fig. (TUl Both models include a dust scattered 
component corresponding to a = OA; g = 0.7 but model A (dashed fine) includes a fiat 
component representing the uncertainty in the baseline and model B includes an extragalactic 
component extincted by the foreground dust. The fit is much better in the case of model 
A with a reduced x^ fit of 0.41 (with 862 degrees of freedom) as opposed to the 0.63 in the 
case of model B and we can formally rule out the presence of any extragalactic light at all 
on the basis of the x^ values. However, our model only grossly represents the true situation 
as, for instance, we have only approximated the effects of multiple scattering and have not 
considered the effects of clumping at all. These are the first observations of the diffuse 
background at such high spatial resolution and further work is required to fully understand 
its distribution. 

The x^ values are lower than one would expect given the preliminary nature of our 
model and are most likely indicative of overlarge error bars. As discussed above, the errors 
are dominated by the 5% uncertainty in the fiat field esti mated through repeate d observations 



of point sources at different positions on the detector (iMorrissey et al.l 120071 ) and may not 
be reflective of the uncertainty in the diffuse background at the larger spatial scales of our 
study. We are currently studying this with many more observations of other locations to 
deconvolve the instrumental effects from true sky variations. 



4. CONCLUSIONS 

We have obtained the highest spatial resolution images of the diffuse UV background 
to date with an effective spatial resolutio n of about 2^ T hese observations are in a region of 



moderate optical depth with r > 1 where ISandagd (119761 ) observed, and correctly identifled 



starlight scattered from interstellar dust. We have obtained G^L£^X observations in both the 
FUV (1350 - 1750 A) and NUV (1750 - 2850 A) bands. After subtraction of the foreground 
airglow and zodiacal light, we were left with about 500 - 800 photons cm~^ sr^^ s~^ A~^ 
in both the FUV and NUV bands. The FUV/NUV ratio increased with increasing FUV 
emission suggesting the presence of H2 fluorescent emission in the FUV ranging up to an 
integrated emission of about 250 photons cm~^ sr~^ s~^ A~^. Our values are consistent 



with those of iMartin et aLl (jl990l ) in their scan of this region; however, we observe a spatial 



variabihty that was not possible with their observations. 



When we link our data with the observations of G251.2+73.3 by iHaikala et al.l (Il995l ). 
we find that the scattered UV light increases linearly with the IR emission for low optical 
depths but saturates at optical depths near unity. This is as expected given that the thermal 
IR emission is from the entire volume of the cloud because of the low optical depth in the 
IR. We have used the same models as in our earlier studies of the diffuse background and 
have found that a dust scattered component with a = 0.4; g = 0.7 is consistent with the 
data; ie., the dust is strongly forward scattering with a moderate albedo. However, the data 
show an intrinsic scatter much greater than can be attributed to photon noise alone which 
must reflect structure at a spatial scale of at least 2' possibly due to variations in the ISRF 
and in the distribution of the interstellar dust. 

We are now extending our analysis to a much larger body of GALEX observations, both 
our own and archival data. Such an investigation will help resolve some of the uncertainties 
in this work such as the contribution of airglow and the zodiacal light. Previous observations 
of the scattering by interstellar dust were on much larger spatial scales and indicated a local 
origin to much of the background; ie., both nearby hot stars and interstellar dust were 
required. GALEX observations will allow us to probe the diffuse background at much higher 
spatial resolutions and thus to investigate the small scale structure of the ISM. 

This research is based on data from the NASA's GALEX GI program, GALEX GIl- 
005007- J092810p702308. GALEX whch is operated for NASA by the California Institute of 
Technology under NASA contract NAS5-98034. We have also made use of NASA's Astro- 
physics Data System and the SIMBAD database operated at CDS, Strasbourg, France. We 
acknowledge the use of NASA's Sky View facility ( |http://skyview. gsfc.nasa.gov ) located at 
NASA Goddard Space Flight Center. 

NVS is supported by a DST Young Scientist award. Support for RCH and LB was 
provided by NASA GALEX grants NNG05GF19G and NNG06GF53G to the Johns Hopkins 
University. 

We most gratefully thank Patrick Morrissey for extensive discussions and clarifications 
beyond the call of duty. 
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Table 1. OBSERVATION LOG 





FUV 


NUV 


Wavelength range 


1350 - 1750 A 1750 - 2850 A 


^eff 


1540 A 


2320 A 


FOV 


1.27° 


1.24° 


Image resolution 


4.2" 


5.3" 


Exposure Time (s) 


14,821 


35,210 


Number of Visits 


10 


22 


RA 




09 28 07 


Dec 




70 21 26 


1 




142.3 


b 




38.2 
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Table 2. FUV AIRGLOW DETAILS 



Visit FUV Minimum 



AGc' 



AGr^ 



'V 



(photons cm ^ sr " s 



Total AG 

-1 .-1 l-i) 



1 


915 


45 


66 


111 


2 


921 


45 


70 


115 


3 


918 


50 


49 


99 


4 


917 


50 


44 


94 


5 


930 


55 


50 


105 


6 


912 


55 


65 


120 


7 


861 





45 


45 


8 


869 


5 


51 


56 


9 


862 


10 


46 


56 


10 


866 


5 


32 


37 



Average Integrated Airglow 



^Relative to Visit 7 on Jan. 3, 2006. 
^Area under the curve above AG^- 



-13- 



Table 3. CONTRIBUTORS TO DIFFUSE EMISSION 

Component FUV NUV 

(photons cni"^ sr"-*^ s~^ A~^) 



Dark Count 


5 


5 


Airglow 


85 


120 


Zodiacal Light 


- 


445 


H2 Fluorescence 


0-250 


- 



Total 90 - 340 570 
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Fig. 1. — Coadded images of FUV (top) and NUV (bottom) intensities observed by GALEX. 
Note the very few and faint stars in both images. 
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Fig. 2. — Total count rate (TEC) converted into photons cm ^ sr ^ s ^ A ^ in the NUV 
channel as a function of local time for each of the visits. 
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Fig. 3. — Effective TEC level in each visit versus solar flux at the Earth. 
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Fig. 4. — G ALEX FUV (top) and NUV (bottom) diffuse images (with point sources removed) 
at a spatial resolution of 2'. The scale is in units of photons cm~^ sr~^ s~^ A~^. 
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Fig. 5. — FUV surface brightness with NUV surface brightness. Airglow has been subtracted 
from both channels with an additional zodiacal light subtraction from the NUV band. Rep- 
resentative error bars equal to about 5% of the data values are shown. The linear correlation 
coefficient is 0.4 but, as noted in the text, we do not find a linear correlation. 
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Fig. 6. — Observed ratio between the FUV and NUV bands with representative error 
bars. The Unear correlation coefficient is now 0.7. 
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Fig. 7. — The predicted level of H2 emission in photons cm ^ sr ^ s ^ A ^of the field. IRAS 
100 /im intensity contours in MJy sr^^ are overplotted and marked. 
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Fig. 8. — The correlation between the IRAS 100 /im intensity in MJy sr ^ and the H2 



emission in photons cm ^ sr 
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is shown. With a correlation coefficient of -0.1, there 



is clearly no correlation between the two. 
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Fig. 9. — IR-NUV correlation. The solid line is the observed slope of iHaikala et al.l (119951 ) 
from Mf/^r observations of the cloud G251.2+73 for the FUV (1400 - I8OOA) region. 
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Fig. 10. — Observed diffuse light with airglow and zodiacal light subtracted ('+') is plotted 
against the optical depth in the region. Overplotted lines are (a) model prediction of dust 
scattered star light (solid line), (b) dust scattered light plus a flat component representing 
an undersubtraction of airglow and zodiacal light (Model A - dashed line) (c) dust scattered 
light plus an extragalactic component (Model B - dotted line). A representative error bar 
for the NUV data is shown. 



